What role pressures play to determine the final end-state of gravitational collapse? 
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We examine here in what way the pressures affect the finaf fate of a continual gravitational 
collapse. It is shown that the presence of a non-vanishing pressure gradient in the collapsing cloud 
can determine directly the epoch of formation of trapped surfaces and the apparent horizon, thus 
changing the causal structure in the vicinity of singularity. 
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Can non-zero pressures within a collapsing matter cloud avoid the naked singularity forming as end-state of a continual 
gravitational collapse? There have been speculations over past decades that when the effects of pressure within a 
collapsing cloud are carefully taken into account, may be only a black hole will form as the end product of collapse 
(see e.g. [1], for an early mention to role of pressure in collapse, and [2] and [3] for further details and references). 
It is clear that in the final stages of an endless collapse pressures will certainly become quite important. Though 
many gravitational collapse models are known currently with non- vanishing pressures present, which end up in the 
formation of a black hole or a naked singularity (see e.g. [4], [5] for a recent discussion), the actual role the pressures 
play to determine the end state of a continual collapse is however not clearly understood. 

We analyze this issue here in some detail in order to bring out the role pressures can play towards determining 
the end state of collapse. A spherically symmetric collapse model is considered where the initial density and radial 
pressure distributions are chosen to be homogeneous, but the tangential pressure is allowed to have a non-vanishing 
gradient, in order to see in a transparent manner the effects it can cause on the evolution of the cloud, and eventually 
towards determining the final state of collapse. It turns out that this by itself can cause inhomogeneities in the density 
, evolution as the collapse develops, thus deforming the trapped surface formation within the cloud. Hence it is seen 
\D 1 that the presence of a non-vanishing pressure can create either of the naked singularity or a black hole as the final 
state for the cloud. Further, the pressure diverges along various families of non-spacelike curves terminating into the 
naked singularity. This is relevant because if a naked singularity developed but if the pressures remained finite in the 
limit of approach to the same, this may not be regarded as a physically interesting situation. 
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The spherically symmetric metric in a general form can be written as, 
CT' 

U~. ds 2 = ~e 2 ^dt 2 + e 2 ^dr 2 + R 2 (t,r)dn 2 (1) 

where dfl 2 is the line element on a two-sphere. Choosing the frame to be comoving, the stress-energy tensor for a 
• i-H , general (type I) matter field is given in a diagonal form [6] , 

X 

Tt = -p; T;= Vr - T e e =Tt= Pe (2) 



The quantities p, p r and p$ are the density, radial and tangential pressures respectively. We take the matter field to 
satisfy the weak energy condition, that is, the energy density as measured by any local observer be non-negative, and 
for any timelike vector V\ we have, 

T lk V l V k > (3) 

which amounts to, 

p> 0; p + p r > 0; p + pe > (4) 

The initial data consists of values of three metric functions and the density and pressures at the initial time t — in 
terms of the six arbitrary functions of the radial coordinate, v{ti,r) — v (r),ip(ti,r) = ip (r), R(ti,r) — r,p(ti,r) — 
Po(r),p r (ti,r) = Prv{r),Pe{ti,r) — pg , where, using the scaling freedom for the radial co-ordinate r we have chosen 
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R(ti,r) = r at the initial epoch. The dynamic evolution of the initial data is then determined by the Einstein 
equations, which for the metric (|l|) become (8irG = c = 1), 

(5) 



r R 2 R" " R 2 R 

y , = 2(Pe-Pr) R' _ P' r , v 

P+Pr R P+Pr 

■ , , G • H , , 

-2i?' + i?'- + R— = (7) 
G H 

G - H = 1 - | (8) 

Here _F = _F(t, r) is an arbitrary function, and in spherically symmetric spacetimes it has the interpretation of the mass 
function, with F > 0. In order to preserve the regularity of the initial data, F(U, 0) = 0, i.e. the mass function should 
vanish at the center of the cloud. The functions G and H are defined as G(t, r) = e~ 2 ^(R') 2 and H (t, r) = e~ 2v (R) 2 . 

All the initial data above are not mutually independent, as from equation (^J) we find that the function voir) is 
determined in terms of rest of the initial data. Also, by rescaling of the radial coordinate r we have reduced the number 
of independent initial data functions to four. We then have a total of five field equations with seven unknowns, p, p r , 
pg, tf>, v, R, and F, giving us the freedom of choice of two free functions. Selection of these functions, subject to the 
given initial data and weak energy condition, determines the matter distribution and metric of the space-time and 
thus leads to a particular collapse evolution of the initial data. 

Consider, for the sake of clarity, the following choice of the allowed free functions, F(t, r) and v{t, r) (see also [7]), 

F(t,r) = lr 3 -±R 3 (9) 

and, 

u(t,r) = c{t)+va(R) (10) 

Also we write, 

R(t,r) = rv(t,r) (11) 

where, 

v(ti,r) = 1; »(t a (r),r) =0; v< (12) 



Using equation ^ in equation (||), we get, 



1 77 ~ I? Pr = l (13) 



v 2 (v + rv') 



It is clear that p(ti,r) = po{r) = 1, i.e. at the initial epoch the density is homogeneous, and as v — * 0, p — ► oo. 
Hence, at the singularity p becomes infinite. Also, we note that the radial pressure remains constant throughout the 
collapse. However, this need not be the case for the tangential pressure, which depends via the Einstein equation (8) 
on the choice of the function v. Then using equation(|13|) we have 

Mr) - f ( P ~^) dr (14) 



Iq \ r 

Assuming that pressure gradients vanish at the center of the cloud, we can take the form of va(r) as, 

Mr)=r 2 g(r) (15) 
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where, g(r) is a suitably differentiable function of r. In that case we can write peo in the following form, 

peo = 1 + r 2 pe2 + r 3 pez H (16) 

where, pg n is proportional to the n th derivative of the initial tangential pressure at the center. We note that choosing 
the form of as given by (|l5|) automatically fixes pei = 0. Also using equation (|l0|) in equation^), we get, 

G(t,r) = b{r)e 2v °W (17) 

Here b(r) is another arbitrary function of r. In corresponding dust models, we can write, b(r) = 1 + f(r), where /(r) 
is the velocity distribution function of the collapsing shells. In the marginally bound case, f(r) = 0. We choose the 
similar analog here and henceforth consider b(r) = 1. 

The reason for the choice such as above for the mass function, and the function b(r), is to bring out the role of 
pressure towards determining the final fate of collapse in a transparent manner. For example, in the present situation, if 
the tangential pressure were vanishing identically, the density evolution would be necessarily homogeneous throughout, 
and the collapse will necessarily end in a back hole, just as the Oppenheimer-Snyder homogeneous dust cloud collapse. 
On the other hand, a non- vanishing pressure gradient gives rise to either one of a naked singularity or a black hole as 
we show below. One can match the cloud to an exterior by introducing an in between shell, wherein p r tends to zero 
at the outer boundary of the shell. 

Using equation(|i~7|) in equation^), we get, 



VRR = -a{t)e 2 " o{R) J R 3 h(R) + |r 3 - \r 3 (18) 
V o o 

Here a(t) is a function of time. By a suitable scaling of the time coordinate, we can always take a(t) = 1. The negative 
sign is due to the fact that R < which is the collapse condition. The function h(R) — h(rv) is defined as, 

p 2u (rv) _ i 

k (™) = r 2 v 2 ( 19 ) 

Substituting the value of the above equation can be written as, 

2 

h(rv) = po2 + -^rvp e3 H (20) 



Now simplifying equation(|l8|), we get 



ni , = - e 2 Mrv)J v 3 I h ( rv j (21) 



Integrating the above equation, we have, 



/«•■<>=/ (22) 

1] 



v ije*"° [v 3 (h(rv) - ±) + 



We note that the coordinate r is treated as a constant in the above equation. Expanding t(v, r) around the center, 
we get, 

t(v,r) = t{v,0)+rX(v) + O(r 2 ) (23) 

where the function X{v) is given by, 

X{v) = -\tdv . V4 ^ P93 (24) 



Thus the time taken for the central shell at r = to reach the singularity is given by. 

f 1 ^/vdv 



UU)= I ^ d \. (25) 
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From the above equation it is clear that for t s (0) to be defined, pe2 > 1/3. Also, the time taken for the other shells 
to reach the singularity can be given as, 

t s {r) =t s (0) + rX(0) + O{r 2 ) (26) 
To see the dynamic evolution of pg(t, r), we put equation (|l3|) in equation(^) and simplify to get, 



pe(r, v) = 1 + - (peo(rv) - 1) (p(r, v) + 1) 



(27) 



It is evident that if at the initial epoch pgo = 1, then pg(r,v) = 1 throughout the collapse, and the evolution will be 
exactly like dust models. In the case otherwise, we find the limiting value of pgo at the central shell as, 



lim lim pg(r, v) = lim lim 1 H 

v— »Q r— >0 v — >0 t — >0 \ v 



(28) 



It is easily seen that near the central singularity there exist ingoing families of timelike curves of the form, 

t S0 -t = kr l3 -k>Q,(3>2 (29) 

along which the tangential pressure necessarily diverges as we approach the point (t So , 0) on the (t, r) plane. Hence 
there always exist some timelike paths along which both the tangential pressure and the density diverge in the limit 
of approach to the central singularity. Thus, as opposed to certain cases, where the pressures are necessarily bounded 
at the naked singularity (e.g. in the case of [1]), which is somewhat artificial situation, we deal here with a physically 
more realistic scenario where pressures diverge at the singularity and then it is to be seen if a naked singularity is 
allowed in such a situation. 

In order to decide the final fate of collapse in terms of either a black hole or a naked singularity, we need to study 
the behaviour of the apparent horizon, and to examine if there are any families of outgoing nonspacclikc trajectories, 
which terminate in the past at the singularity. The apparent horizon within the collapsing cloud is given by R/F = 1, 
which gives the boundary of the trapped surface region of the space-time. If the neighborhood of the center gets 
trapped earlier than the singularity, then it will be covered and a black hole will be the final state of the collapse. In 
the case otherwise, the singularity can be naked with non-spacelike future directed trajectories escaping from it to 
outside observers. 

To consider the possibility of existence of such families, and to examine the nature of the central singularity occurring 
at R = 0, r = in the present case, let us consider the outgoing radial null geodesies equation, 



dt 
dr 



(30) 



The central singularity occurs at v = 0,r = 0, which corresponds to R = 0,r — 0. Therefore, if we have any future 
directed null geodesies terminating in the past at the singularity, we must have R — > as t — > t s (Q) along the same. 
Now writing the geodesic equation equation j30|) in terms of the variables (u = ri , R), we have, 



dR 3 _ 2 , 
du 5 



1 + §e+-» 



(31) 



Using equation(0) and considering R < 0, we get, 




(32) 



If the radial null geodesies terminate at the singularity in the past with a definite tangent, then at the singularity 
the tangent to the geodesic 4^ > 0, in the (u, R) plane, with a finite value. In the present case, all singularities for 

-co. Therefore, only the singularity at the central 



r > are covered since ^ 



oo in that case, and hence 4^ 



shell could be naked. Now from equation (pl|) we get for r — > along a constant v line, 



X(v) 




P92 



(33) 
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Let us define the tangent to the null geodesies from the singularity as 

R dR 



xq = lim lim 

t^t B r^o u du 



(34) 

t— >*.:) — >0 



Using equation (B2), we get, 



4 = -^X(O) (35) 



Now, if pg3 < 0, then Xq > and hence we would have radially outgoing null geodesic coming out from the singularity, 
and the singularity will be naked. While if pg 3 > 0, we will get a black hole solution. If pg 3 = 0, then we will have to 
go to the higher order terms and do the same analysis. In the (t, r) plane, the equation for the radial null geodesic 
coming out from the singularity is, 

t-t s (0) =x Q ri (36) 

Also, we see that in case of a naked singularity, the singularity curve at the center, equation (p3[), is an increasing 
function of r, as in that case X(v) > 0. On the other hand, a black-hole solution gives a decreasing or constant 
curve for shells with increasing r. It is relevant to note here that X(0) > implies v' > 0, and so from R' = v + rv' 
we see that there are no shell-crossings, at least in a neighbourhood of the central singularity. We thus see how a 
non- vanishing pressure causes a naked singularity as the end state of collapse. 

It is also interesting to note how the non- vanishing pressure gradient affects the spacetime shear. From equation([l3]), 
we get for small values of r along constant v line, 



p(t, r) = , - 1 (37) 

3/^3 



If P83 < 0, then X(v) > 0, i.e. v' > 0. Thus from the above equation it is evident that p(t,r) is a decreasing 
function of r at any given time t. But it is known [8] that in the case of a vanishing shear, for the mass function we 
have considered here in equation(||), p = p(t), i.e. the density is necessarily homogeneous throughout the collapse. 
Therefore, we conclude that in the case of a naked singularity developing as collapse end state, the collapsing cloud 
has non-zero shear. This shear may play the role of deforming the apparent horizon, thus exposing the singularity. 
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